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НЕЛИНЕЙНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ 
ГИДРОДИНАМИЧЕСКОЙ СМАЗКИ СЛОЖНОНАГРУЖЕННОГО 
СОСТАВНОГО КОНИЧЕСКОГО ПОДШИПНИКА 

С ПОРИСТЫМ СЛОЕМ НА ЕГО РАБОЧЕЙ ПОВЕРХНОСТИ 


В работе на основе нестационарных нелинейных уравнений Навье-Стокса для слу- 
чая «тонкого слоя» и уравнения Дарси с использованием принципа конструктивной 
суперпозиции разработан метод гидродинамического расчета составного кониче- 
ского подшипника. Дана оценка влияния угла конусности и протяженности пори- 
стой составляющей на основные рабочие характеристики составного подшипника. 
Ключевые слова: математическая модель, гидродинамическая смазка, составной 
подшипник, пористый слой. 


Введение. Как известно [1], конические подшипники жидкостного трения 
широко применяются как в тяжелонагруженных узлах машин и механиз- 
мов, так и в качестве конструктивного элемента использования опорных 
узлов высокоскоростных роторных машин. Применение конических ра- 
диально-упорных подшипников в узлах высокоскоростных опорных машин 
позволяет уменьшить габаритные размеры опор, снизить энергетические 
потери, улучшить динамические характеристики опорных узлов [1,2]. 
Поэтому разработка метода гидродинамического расчета конических под- 
шипников представляет теоретический и практический интерес. Анализ 
публикаций в данной области показывает, что в работах, посвященных 
расчету конических подшипников, рассматриваются либо сплошные кони- 
ческие подшипники [1], либо содержащие пористый слой на всей рабочей 
поверхности [3]. В данной работе приводится метод гидродинамического 
расчета составного конического подшипника, работающего в нестационар- 
ном режиме. 

Постановка задачи. Рассматривается неустановившееся движение вяз- 
кой несжимаемой жидкости в зазоре составного конического подшипника 
конечной длины. Вкладыш, представляющий собой совокупность сплошных 
и пористых втулок, запрессованных в непроницаемый конический корпус, 
предполагается неподвижным, а конический шип вращается с постоянной 


угловой скоростью @, на которую накладываются возмущения @ (Г). Так- 


же предполагается, что смазка подается в осевом направлении и в двух 
торцевых сечениях давление задано. Поместим начало цилиндрической си- 


стемы г,0,2 в левом конце на оси подшипника (рис. 1), ось 2 направим по 


оси подшипника. В выбранной системе координат уравнения контуров ©, & 
и а записываются в виде: 


г=7+2120; Г=И+ 210; ГЕК 2190 + ес0$9, (1) 
где '- радиус шипа; к- радиус вкладыша в сечении 2 =0; 
б = ( м ") с05й - толщина пористого слоя; е -— эксцентриситет; 


@- угол конусности подшипника. 
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Рис. 1. Схематическое изображение составного конического подшипника 


Основные уравнения и граничные условия. Будем исходить из нели- 
нейных нестационарных уравнений Навье-Стокса для «тонкого слоя» [3], 
уравнения неразрывности и уравнения Дарси: 
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Система уравнений (2) решается при следующих предельных усло- 
вИЯХ: 


- прилипание смазки к поверхности шипа и подшипника; 











: (2) 
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- на внутренней поверхности вкладыша [при 2 [= || нор- 
мальная составляющая скорости определяется законом Дарси; 


- при переходе через пористую границу | Ея, | давление 


меняется непрерывно; 

- на внешней поверхности пористого слоя нормальная состав- 
ляющая скорости равна нулю; 

- в двух сечениях (в начальном и конечном) давление задано. 


Приведем начальные условия: 
= (0) _ д (0) _ д (0) _ д (9) = (0) 
Р|.0=Р 50, [1-070, о 20:11:02: РР - (3) 


Здесь р ® 0 ®у®,Р ©) - решение соответствующей стационарной за- 
дачи. 





Для решения поставленной задачи используется прием конструк- 
тивной суперпозиции, заключающийся в следующем: сначала последова- 
тельно решаются две задачи, связанные с разработкой метода расчета со- 
ставных цилиндрических подшипников с конструктивными параметрами, 
соответствующими сечениям 2 = Ои 2 =[ составного конического под- 
шипника. В последующем определяется линейная суперпозиция этих реше- 
ний. 

Задача 1. Разработка метода расчета составного цилиндрического 
подшипника с радиусом шипа и, и вкладыша и и с толщиной пористой со- 


ГОО) 


ставляющей втулки б. Определение поля скоростей и.’0,’,0’’ и давле- 


ний соответственно в смазочном слое р“ и в пористом Р®. 

Задача 2. Расчет составного цилиндрического подшипника с ради- 
усом шипа х+ (за и вкладыша + Иеа . Определение поля скоростей 
у и полей давлений р” и Р®. С использованием решений за- 


дачи 1 и задачи 2 приближенное решение рассматриваемой задачи пред- 
ставляется в виде линейной суперпозиции этих решений, т.е. [4]: 


* * * * 
р= р 1% а о к -в ), 
9, =009 "+0 2 (1-0 °), 0. =009 "+00)(1-0°), 
Р=Р“9`+Р®(1-9°),8° [01]. (4) 


Найдем решение задачи 1. Перейдем в смазочном слое к безраз- 
мерным переменным по формулам: 


и ОО Е 
= ди’, 9, = 07,0. =0КИ 


Фо 
›Р -РР, 
И И 


2700” 


2 
* и * 
се. ЕЕ ==. (5) 


6 2 
В пористом слое переход к безразмерным переменным осуществим 
по формулам: 
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г * РС 1 * РС 
На: Вр РОУ Ор: ао 
Используя (5) и (6), систему уравнений (2) и граничные условия (3) 
В сЛучае рассматриваемого цилиндрического подшипника с пористым сло- 
ем на рабочей поверхности, можно записать решение задачи в виде (с 


7? 


б 
точностью до членов О — ‚О — ): 
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9” 08 92 91 9ё Е 9 00 92 
д9Р 
а = м на (1) ко = 0, У -0 = 0, р= Р ОИ 
96-2 ий 
И. 0 1510 е 
ем = -1 518 - Е 
636 7=1-1 с0$8 И 0 1 б 
0 
р - 7=1-1 с0$8 = Те Ри. и 7=1-1 с058 = 0, 
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Ро Ир р 192 Гр, 
Р® т: 2 
| =0, №: бпри2 А 7 ее. 
а Ей 1 Г 


(1) = и) (1) < ›,(0) (1) = д (0) ао (1) < р) 
р |6=Р ‚и |= и 0 |0=% У нш=и,Р аНЫР . (8) 


0 0 0 0 0 к о 
Здесь р. и у и Р® - решение соответствующей безразмерной 
стационарной задачи. 
В системе уравнений (7) осредним по толщине смазочного слоя 
члены, обусловленные силой инерции, а также осредним уравнение Дарси 
по толщине пористого слоя. Введем следующие обозначения: 


1-1 0998 4 (1) (о ( (0 





+ Ве и +0 + м а, 
1-1 с080 р 91 ди 90 92 
1-1 с0$0 
МО = Ре. 3 | тв и + Ве и? ры и +0 оби и + в и” Чи 
1-1 с080 т 91 ди 99 92 . 
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В (1) Г 2 ра) (1) 2ра) 
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Точное автомодельное решение задачи (7) - (8) с учетом (9) будем 
искать в виде: 


(0 0) 


их = Е.В), ИФ = и (Е, В (8), ро. и”, 0), 

















в" Я 6,0, 2 мов да = р шв 
РО = 490,6) м На + 49) —- 1-5 —. 1 5108 .(10) 
Подставив (10) в (7) и (8) с и ее будем иметь: 
ну © 
ЭР. = & (0), ЕЕ а тт: 011) М1? (8 ), 
9 \ 9 и” в“ я“ *(1) 
о В В. 
р” 94% 92 40 92 р 
на =." 90 2 + А (В+ В,)+ — (11) 
(1) 
ар’ 97 _ А“ 
= 156) =0, и 9,г)= - 
ое и) ое (= @,0= 
00,0) = 0, #0, = 0, а, = 1- . 
() 
-@) оо м ор 
| (1,1) = -1+ А ВАЙ (1,1) = 0, р = Р ри „= 0, 
(7 т МЕ 
Аб (0,5) = АЗ(2л 1); ро (0,24) = РО(21 ‚24. — (2) 
Здесь 
аа бе бое 
бт в На. ИЕ а -- 4; 
В 
р р 
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Решение задачи (11) - (12) находится непосредственным интегрированием. 
В результате будем иметь: 
(1) 
9 _ - 2 2_ я 9% 
6 2.05 
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(1) 














о ми 
И Ш в} 
9 о 9 
р® = с - 6 (6) @ + а2+6- М!’ ,2)4, 
ой) #9) о 
м8 6 =-—2- т: 148, 
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Е д 
т 9ё р) 1? 9 
т 1-ю ВУ 190 в ой 
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_ Ва В5-Вь, в о 0 о 
А (те 8" А’, а = РВ +12М0 4, А = сои. 


Явный вид функции А? ‚Ё) моб ‚Г) при определении основ- 
ных рабочих характеристик нам не понадобится. Соблюдая периодичность 
гидродинамического давления по 60, для определения &,(!) приходим к 


следующему уравнению: 
21 21 


а ; н #- М%6@,0а = 0. (14) 
0 0 


Таким образом, задача 1 полностью решена. Задача 2 решается по 
той же схеме. Переход к безразмерным переменным осуществляется по 


формулам (5) с заменой в рассматриваемом случае х, на и, + Пэй ик на 
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(2) -2) 20 
п+ аа. ДлЯУу и т 7 и ре Р@) получим выражения, аналогич- 


ные (13) и (14). 

Перейдем к определению воздействия смазки на подшипник. Для 
проекций главного вектора и главного момента в случае задачи 1 будем 
иметь следующие выражения: 








272 
ВОР рзв 48 4, 
д 00 
у уро ен т 
р —+——— 0442, Мри® а 
2 ры 92 |2 9Ё й Е ря п п о 
В случае задачи 2: 
3 721 
ВО ВВМ рэ 4 4, 
б 00 
у ро (2) ^ (2) 
МР = в (+ Ива (+ аа и а Чт 46, 
т 
$=0 
2. 2) 





М№йри@ (16) 


2 
+ Поа к+ Поа 
Используя метод суперпозиции для основных рабочих характери- 
стик, получим выражения: 


_ рад * (2) * 
В, = 88 "+ В 8"), 
= (Пр * (2) * ь 
М. = М9 "+ М9(1-6°), 8° [0,1 
Из найденных аналитических выражений ДЛЯ ОСНОВНЫХ рабочих ха- 


рактеристик составного конического подшипника следует, что эти характе- 
ристики существенно зависят от следующих параметров: 


мо НР 


скорость на рабочей пористой поверхности составной цилиндрической 


(17) 


(2) 
‚№ А’, соответственно характеризующих безразмерную 


втулки радиусом И и радиусом И + Пой ; 0/0 - характеризующего 


относительное возмущение, накладываемое на угловую скорость вращения 


1 1 К НИ 
вала; м! МО — обусловленных силой инерции смазочной жидкости; 
=) 2 А 
В” - обусловленного нестационарностью уравнения Дарси; РЕ 


характеризующего протяженность пористой компоненты составного кони- 
ческого подшипника; ( - характеризующего конусность подшипника. 
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Обсуждение результатов. Результаты численного анализа (рис.2) пока- 
зывают следующее: 

1. Безразмерная несущая способность и безразмерный момент тре- 
ния с увеличением значения безразмерного времени стремятся к значени- 
ям, соответствующим стационарному режиму. 

2. Наиболее оптимальным (по несущей способности и моменту тре- 
ния) является значение угла конусности, близкое к 6°, и относительная 
протяженность пористой составляющей, равная 1/6. 
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Рис. 2. Зависимость безразмерной несущей способности (а) и безразмерного момен- 
та трения (6) от безразмерного времени при различных значениях угла конусности 
и относительной протяженности пористой составляющей конического подшипника: 


*() =(1) 


МА = 0,2: В 


*(1) 


=:10: 0/0: = Тао Ро“ = бе == 


1- -а=6,, Е Ш: э И 
Па 4: © п 6 
И Е 
п 2 п 6 

- СПЛОШНОЙ ПОДШИПНИК, - - - -- составной подшипник 
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Заключение. Предложена нелинейная математическая модель гидроди- 
намического расчета составного конического подшипника. Установлены оп- 
тимальные (по несущей способности и моменту трения) значения угла ко- 
нусности и протяженности пористой составляющей составного конического 
подшипника. 

Обозначения: р — гидродинамическое давление в смазочном слое; 
Р - гидродинамическое давление в пористом слое; е — эксцентриситет; 
7’, — радиус шипа; и — радиус вкладыша (в сечении 2=0); 7 - и — толщина 
пористого слоя; и,ё{,2 - цилиндрические координаты; Р -— плотность 


1 К 
смазки; К°” — проекция главного вектора воздействия смазки на подшип- 
2. Я 
ник; М® - момент силы трения в случае задачи 1; А”, М - то же в 


случае задачи 2; 9,0,,0, — компоненты вектора скорости; К — проницае- 
мость; { -— динамический коэффициент вязкости. 
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п {15 гезеагсИ а те оЧ оЁ Не Пудгодупатюс сасшаЧоп оЁ {Пе сотрозке 
сопс Беаипд 15 ммогкеЯ оц оп Ба$!5 оЁ {Ле поп-5аЧопагу едицаНоп$ оЁ Мамег- 
ЗюкКеб Гог {Пе сазе оЁ {Пе «т |ауег» апа Не Багзу едиаНоп м Ше изе оЁ 
{Не те!поЧ оЁ {пе соп&гисНопа! зирегро$оп. Т№М5 гезеагсН ргезетс5 ап ар- 
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